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ABSTRACT 

Previous surveys in a few metal-poor globular clusters (GCs) showed that the determination of 
abundances for Li and proton-capture elements offers a key tool to address the intracluster pollution 
scenario. In this Letter, we present Na, O, and Li abundances in a large sample of dwarf stars in the 
metal- rich GC 47 Tucanae. We found a clear Na— O anticorrelation, in good agreement with what 
obtained for giant members by Carretta et al. While lithium and oxygen abundances appear to be 
positively correlated with each other, there is a large scatter, well exceeding observational errors, and 
no anticorrelation with sodium. These findings suggest that Li depletion, due to mechanisms internal 
to the stars (which are cooler and more metal-rich than those on the Spite plateau) combines with 
the usual pollution scenario, responsible for the Na— O anticorrelation. 

Subject headings: globular clusters: individual (NGC 104) - Stars: abundances - stars: individual 
(Population II) 



1. INTRODUCTION 

The traditional paradigm of globular clusters (GCs) as 
the textbook example of a simple stellar population is 
now outdated. In the last several years, a wealth of ob- 
servational studies (both photometric and spectroscopic) 
have been carried out, revealing the quite complex nature 
of this old Galactic population. 

With the famous exception of Cen (Freeman & 
Rodgers 1975), and of the more recently scrutinized M 54 
(Carretta et al. 2010), M22 (Marino et al. 2009), and 
Terzan 5 (Ferraro et al. 2009), GCs show homogeneous 
compositions in the iron peak and heavier a elements 
(e.g., Ca and Ti). Abundance variations in the lighter 
elements, namely Li, C, N, O, Na, Mg and Al, have been 
instead recognized since the 1970's (e.g. Cohen 1978; see 
Gratton, Sneden & Carretta 2004 for a recent review). 
These peculiarities of GC stars require that some ma- 
terial must have been processed through the complete 
CNO cycle in hot H burning, through proton-capture re- 
actions (Denisenkov & Denisenkova 1990): element pairs 
C and N, O and Na, and Mg and Al are anticorrelated, 
the abundances of C, O and Mg being depleted while 
those of N, Na and Al being enhanced. 

A previous generation of stars, which synthesized 
proton-capture elements in their interiors, is now com- 
monly accepted as responsible for the above-mentioned 
chemical signatures. Regardless of the debated nature of 
these element polluters (asymptotic giant branch (AGB) 
stars experiencing hot bottom burning, e.g., Ventura & 
D'Antona 2009, or fast rotating massive stars (FRMS), 
e.g., Decressin et al. 2007), two fundamental observa- 
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tional facts deserve to be emphasized. (1) All the GCs 
surveyed so far show the Na— O anticorrelation (Carretta 
et al. 2009b): this indicates the presence of at least two 
populations (neither coeval nor chemically homogeneous) 
within each cluster. (2) The primordial nature of such 
a phenomenon is uniquely indicated by the presence of 
these chemical signatures also among not evolved turnoff 
(TO) or scarcely evolved (subgiant, SGB) members (see, 
e.g., Gratton et al. 2001). 

In this context, Li abundances play a fundamental role. 
In fact, this element can be easily destroyed in stellar 
interiors (starting at Tburn~2.5xl0^ K); and since the 
CNO and NeNa cycles require much higher temperatures, 
it is expected that in the site(s) where these cycles occur 
no Li is left. In particular, the Na-poor, O and Li rich 
stars, that are the first population born in the cluster, 
share the chemical composition of field stars of the same 
metallicity, while the Na-rich (Li- and 0-poor) stars form 
from gas progressively enriched by ejecta (which are rich 
in Na, depleted in O and Li) of the first generation. 

As a consequence, if primordial and processed mate- 
rial are mixed in different proportions, then Li and Na 
(Li and O) are expected to be anticorrelated (correlated) 
with each other. 

Measuring the Li abundances in unevolved stars pro- 
vides a direct indication of the amount of pristine, and 
by difference of polluted (CNO-cycle processed), mate- 
rial present in each star. This makes Li a unique tracer 
of the dilution process which took place in the GC, sup- 
plying also fundamental insights on its early evolutionary 
stages. Another important point is that Li abundances 
can provide strong observational constraints on the ori- 
gin of the polluters: since AGB stars might also pro- 
duce Li (via the Cameron— Fowler mechanism, Cameron 
& Fowler 1971), while massive stars can only destroy it, 
if Li-rich stars are present in GCs, AGB stars would be 
definitely favored with respect to FRMS. 

To date, only three GCs have been surveyed for correla- 
tions between Li and proton-capture elements. Pasquini 
et al. (2005) obtained Li abundances in nine TO mem- 
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bers of NGC 6752: they found a depletfon reaching down 
to ^1 dex below the Spite plateau values, clearly anti- 
correlated with Na abundances. A similar result was 
obtained by Bonifacio et al. (2007) who found a scatter 
in Li abundances much larger than observational errors 
and anticorrelated with Na among (only) four stars in 47 
Tuc. Very recently, Lind et al. (2009), targeting about 
100 main-sequence (MS) and early subgiant branch stars 
in NGC 6397, detected for the first time a significant 
anticorrelation between Li and Na in this GC. This last 
investigation supersedes an older research by Bonifacio 
et al. (2002) based on only a few stars, and highlight- 
ing the importance of large samples of stars for similar 
studies. 

In this Letter we present Na, O, and Li abundance 
determinations in ~90 unevolved TO stars of the GC 
47 Tucanae, providing the largest database of this kind 
available in the literature so far. 

2. SAMPLE AND ANALYSIS 

We retrieved from the ESQ Archive FLAMES-Giraffe 
spectra (MEDUSA configuration) of 109 unevolved mem- 
bers of the metal-rich GC 47 Tuc (Program 081. D- 
0287; PI Z. Shen). The observations were carried out 
in Service mode in 2008 August-September; gratings 
HR15n (R -17,000), HR18 {R -18,400), and HR20A 
{R -16,000) were used to cover Li i (6707 A), O i (7771, 
7773, 7774 A), and Na i (8183, 8194 A) features, respec- 
tively. 

The data reduction was performed with the ESO Gi- 
raffe pipeline^ using Gasgano and following the stan- 
dard procedure: bias subtraction, flat-field correction, 
wavelength calibration, and optimal spectrum extrac- 
tion. A careful sky subtraction has been done using 
the IRAF^ task skytweak in order to take into account 
wavelength shifts and possible differences in emission line 
strengths between stellar and sky spectra. It is notewor- 
thy that the presence of the moon at more than 90% il- 
lumination along with diffuse clouds caused some of the 
observed frames to be severely contaminated and were 
thus discarded (the night of August 17th). Different ex- 
posures were then co-added to improve signal-to-noise ra- 
tios, obtaining typical values of —50-100 per pixel (near 
the Li I region at -6708 A). 

The abundance analysis was carried out using the 
ROSA code (developed by R. Gratton) and Kurucz 
(1993) model atmospheres; for Na and O, we derived 
abundances from the equivalent widths (EWs) , applying 
corrections due to non-LTE effects (Gratton et al. 1999). 
LTE Li abundances were instead obtained through spec- 
tral synthesis. Stellar parameters were derived as follows. 
(1) We obtained initial effective temperatures (Tcff) from 
the de-reddened {E{B - V")=0.04; Harris 1996) B - V 
colors'" and the calibration by Alonso et al. (1996). 
The adopted values were derived from a relation between 
V magnitude and color-based Teg. This reduces errors 




^"http: //www. eso . org/ sci/data-processing/sof tware/pipeline; 

® IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is oprated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

Data are from unpublished photometry by Y. Momany (private 
communication); they are derived and calibrated exactly as the 
other dat sets in Momany et al. (2003). 
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Fig. 1. — Na— O anticorrelation for 47 Tuc dwarf members; error 
bars come from uncertainties on EWs. A dilution model (Prantzos 
& Charbonnel 2006) is also overplotted as a (red) solid line. 

for individual stars. (2) The surface gravities were es- 
timated from the fundamental relation between M, L, 
and Teff (logs = 4.44 + log(M/MQ) - log(L/L0) + 
4 log Toff — 15.0447), by adopting the bolometric correc- 
tions by Alonso et al. along with a distance modulus 
of 13.34 (Harris 1996), while Mto was calculated in the 
same fashion described in Gratton et al. (2010). 

Finally we assumed for all our sample stars microtur- 
bulence ^=1.00 km s^^ (note that none of the elements 
under scrutiny is strongly affected by ^ values). 

The [X/Fc] ratios were calculated assuming as solar 
references log n(Na)Q=6.21 and log n(0)Q=8.79, val- 
ues obtained in a fully consistent fashion from the anal- 
ysis of the solar spectrum by Gratton et al. (2003), and 
[Fe/H]=— 0.76 (Carretta et al. 2009a) as cluster metal- 
licity. 

Both internal (random) and external (systematic) er- 
rors may affect our analysis. A discussion of the first 
kind of uncertainty, which is predominant for investigat- 
ing star-to-star variations, can be found in Carretta et al. 
(2007); here we just mention that errors in EWs are dom- 
inant, while the internal uncertainties on stellar param- 
eters (Toff) can be neglected (note that the distance of a 
star on the color-magnitude diagram from the MS main 
loci represents the major error source in the adopted tem- 
perature values). Systematic errors due to the adopted 
set of model atmospheres and/or to reddening and Tos 
scale cannot be larger than —0.1 dex (see Bonifacio et al. 
2002, 2007 and Gratton et al. 2001); in particular, for 
Li, the combination of three-dimensional and non-LTE 
effects yields negligible corrections (see, e.g., Asplund et 
al. 2003). 

iWr«^ie>^frsH^^Mcpfe'M^itl^ ^'^^ corresponding un- 
cci iaiuLics, aiuiig wlLii niagmiudes and coordinates, are 



available only in electronic version. 

3. RESULTS AND DISCUSSION 

As shown in Figure [U the unevolved cluster stars re- 
veal a very clear Na— O anticorrelation, consistently with 
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what has been obtained for giant members (Carretta et 
al. 2009b). We detected an oxygen variation of about 
'^0.80 dex and a change in Na of ~0.77dex, to be com- 
pared with ~1 and 0.82 (O and Na respectively) for gi- 
ants in Carretta et ah (2009b). The interquartile ranges 
IQR[Na/0] are 0.368±0.067 and 0.470±0.041, respec- 
tively. These two estimates are consistent at Icr level; 
however since the dwarfs are a more homogeneous sam- 
ple (i.e., very similar atmospheric parameters) and the 
spectral features used to measure a given element are not 
the same in dwarfs and giants, the latter being weaker 
and plagued by higher uncertainties, both the spread 
(maximum-minimum) and the IQR result slightly higher 
for giants. It is noteworthy that we found an offset in Na 
abundances between dwarfs and giants of ^^0.2 dex for 47 
Tuc. A similar offset can be obtained also for NGC 6397, 
when comparing Na abundances for dwarfs from Lind et 
al. (2009) and for giants from Carretta et al. (2009b). 
We investigated the nature of this offset, concluding that 
it is due to the adopted set of lines, with the doublets 
5682- 5688 A and 6154-6160 A (used in Carretta et al.) 
yielding always higher abundances than the 8183-8194 A 
features (employed for dwarfs). We, in fact, detected the 
same effect in the Sun (with the same A(Na)'^ 0.2 dex). 
We propose that the way damping wings are treated is 
responsible for the observed difference. The present anal- 
ysis uses values consistent with the Barklem et al. (2000) 
values; however, EWs are extracted using Gaussian pro- 
files, which may result in systematic underestimates of 
the EWs for very strong lines, like the 8183-8194 A ones. 
We then consider abundances from the weaker doublets 
as more reliable, and we offset our Na abundances for 
this effect. 

In general, we conclude that there is a very good agree- 
ment between the two distributions as can be clearly seen 
in FigureO after applying the derived offset in Na abun- 
dances, the cumulative [Na/0] functions look very simi- 
lar, and a Kolmogorov— Smirnov (KS) test indicates that 
the two distributions cannot be distinguished (the prob- 
ability from KS being ^14%). 

Thanks to the wide sample available, we confirm that 
evolutionary effects, acting during the red giant branch 
(RGB) phase (see, e.g., D'Antona & Ventura 2007 for 
M 13), cannot contribute to the Na— O distribution (at 
least for the present cluster), since the extent of anticor- 
relation in both giants and dwarfs is essentially the same. 
Furthermore, following the approach by Carretta et al. 
(2009b) the fractions of primordial (P)* , intermediate (I) 
and extreme (E) stars computed from unevolved mem- 
bers are 34%±5%, 63%±7%, and 3%±1% respectively 
to be compared with 27%±5%, 69%±8%, and 4%±2% 
from giants (see Carretta et al. 2009b). 

On the other hand, there is not a one-to-one correla- 
tion between Li and O abundances (left-hand panel of 
Fig. [3]): even if some trend might be present, there is 
a considerable scatter, much larger than observational 
errors. The figure shows that, while the stars with 
low O values have also low Li content, at higher O 
abundances, Li can assume a large range values, from 
1.54±0.06 to 2.51±0.18, with the most extreme star (# 

* The P stars must have [Na/Fe]<[Na/Fc]min 4- 0.3, then the 
remaining stars are divided in I and E according to [O/Na] larger 
and smaller than —0.9 dex, respectively. 
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Fig. 2. — Cumulative [Na/O] distribution for dwarfs (red solid 
line) and giants (blue dashed line). 

085) having log n{U) = 2.78±0.08 (this very Li-rich 
star, which appears as primordial type according to its O 
abundances, surely deserves special attention and further 
high-resolution observations, in order to investigate this 
extreme behaviour). In Figure |4] we show, as explicative 
example, the spectra of three stars with almost the same 
O abundances; one is star # 085 and the two others have 
intermediate and low Li. 

If we then focus on the Li— Na diagram (middle panel of 
Figure [3]), the scatter is even larger with no evidence for 
an anticorrelation. As further evidence, the linear corre- 
lation coefhcient is 0.39 (77 stars) and —0.02 (84 stars), 
respectively, for the Li— O and Li— Na distributions: in 
the first case, the significance level is >99.5%, while the 
second one has no statistical meaning. 

Our study increases the sample of Li determination in 
this cluster, completing the previous studies by Pasquini 
& Molaro (1997, Li for two stars) and by Bonifacio et 
al. (2007), who proposed a Li— Na anticorrelation on the 
basis of only four members. 

To get more insight into this issue, we can compare it 
with a simple model. It is known that the Na— O an- 
ticorrelation can be well described by a simple dilution 
model^, like the one proposed by Prantzos & Charbonnel 
(2006). On the other hand, Li abundances show a wide 
scatter, with values often well below those expected from 
a simple dilution model (red solid line in three panels of in 
Figure 131). However, we cannot exclude that this model 
might represent the upper envelope of the Li distribu- 
tion. In the case of ^7 Tuc, Li is telling us something 
different from Na and O. Different mechanism(s) have 
to be invoked in order to explain the observed abundance 
pattern. 

The distribution of Li abundances for 47 Tuc hence 
looks very different from that found for the metal-poor 

^ The basic equation of this model is: [X] = log[{l — dil) X 
lOl^o] + dil X lOl^pl], where < dil < 1, [Xq] and [Xp] are the 
logarithmic abundances of the element in the original and processed 
material, respectively. 
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([Fe/H]=-1.99; Carrettaet al. 2009a) GC NGC 6397. In 
that handful of Na-rich, Li-poor stars indicates a 

significant Li— Na anticorrelation (Lind et al. 2009), and 
a simple dilution model very well reproduces the Li and 
Na distributions (see the right-hand panel in Figure [S]) . 

It is important to keep in mind that TO stars in 47 Tuc 
are cooler (Tefi~ 5700-5800 K) than the TO stars of NGC 
6397 (Teff- 6100-6300 K), and they are well below the 
limit usually considered for the Spite plateau. The scat- 
ter found in Li abundances of 47 Tuc stars is reminiscent 
(or even better, is the Population ii analog) of the large 
scatter of Li abundances obtained for solar twins in the 
old open cluster M 67 (Randich et al. 2000) and in gen- 
eral, for old thick disk stars with effective temperatures 
close to '^5800 K (Ryan et al. 2001). We note, however, 
that in this case the Li abundance variation is smaller 
with respect to typical values detected in Population i 
stars, probably because the latter ones are significantly 
more metal- rich. 

4. SUMMARY AND CONCLUDING REMARKS 

We present in this Letter Li, Na, and O abundances for 
a large sample («90) of TO stars in the metal-rich GC 47 
Tuc, providing the largest database of this kind available 
so far. Our main results can be briefiy summarized as 
follows: 

1. We obtained a very clear Na— O anticorrelation, 
confirming the previous findings from giants de- 
rived by Carretta et al. (2009b). This is the first 
time that the Na— O distributions in dwarfs and gi- 
ants can be directly compared for a cluster using 
large samples of comparable sizes for both evolu- 
tionary stages. At least for 47 Tuc, evolutionary 
effects due to the RGB phase can be ruled out as 
contributors to the extent of Na— O anticorrelation, 
which in both cases span the same range (within 
the observational uncertainties). 

2. As expected from stellar nucleosynthesis models in 
conjunction with multiple population scenarios, Li 
abundances should be positively correlated with O 
and anticorrelated with Na. At variance of the 
metal-poor GC NGC 6397, in the case of 47 Tuc 
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the Li content does not show an anticorrelation 
with Na, and only a weak correlation appears with 
O, with a quite scatter distribution from both di- 
agrams. Our result disagrees with the previous 
study by Bonifacio et al. (2007), who found Li— Na 
anticorrelation from a small sample of only four 
stars, and once again emphasizes the crucial role of 
statistics in this kind of analysis. A simple dilution 
model fails to reproduce the Li— Na— O distribu- 
tions for this cluster and advocates the presence of 
some different mechanisms responsible for the ob- 
served Li pattern. 

3. The scatter we find in Li abundances reminds of 
what has been detected, and reported in a large 
body of the literature, in Population I stars of sim- 
ilar parameters {Tcs, \ogg), the most famous case 
being the old open cluster M 67 (see, e.g., Randich 
et al. 2000). 

We are not presently able to conclude if the trend we dis- 
covered in 47 Tuc is peculiar or, on the other hand, other 
GCs share a similar behavior. In fact, to date, only two 
GCs have been investigated from this point of view. In 
this context, we mention that we cannot explain the un- 
likeness in the Li— Na distributions between NGC 6397 
and 47 Tuc; in particular, we cannot discriminate if the 
differences in Teg 's or metallicity can account for such a 
discrepancy. To probe this issue, it is crucial to enlarge 
the sample of simultaneous determinations of Li, Na, O 
in GC dwarf stars, by including other (nearby) clusters 
with different structural parameters (e.g., HB morphol- 
ogy, age, metalhcity). 
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Fig. 3. — Run of Li with O and Na for 47 Tuc (left and middle panels). The (blue) open symbols are stars with Li from Bonifacio et 
al. (2007) and Na, O from Carretta et al. (2004). A dilution model is overplotted as (red) solid lines (see the text). The arrow represents 
the effect on abundances of a change in Toff of 145 K, which is the offset between the two different temperature scales employed. In the 
right-hand panel we show, for comparison, the Li— Na diagram for the GC NGC 6397 from Lind et al. (2009). 




Fig. 4. — Two portions of the spectra, centered on Li (left panel) and O (right panel), are shown for the stars #068, #083, and #085. 
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TABLE 1 
Properties of target stars. 
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Note. — We report our ID in Column 1, R.A. and Decl. in Columns 2 and 3, while magnitudes V and B are given in Columns 4 and 
5. Stellar parameters (Tgg- and logg) along with abundances and their errors are listed from Columns 6 to 13, respectively. (This table 
is available in its entirety in a ma<;hine-readable form in the online journal. A portion is shown here for guidance regarding its form and 
content.) 



